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ABSTRACT 


An- attempt is raade to explain the g Cephei variables as 
non-rotating stars undergoing radial oscillations, on the basis 
of their relatively low observed rotational velocities and the 
period analyses by van Hoof. Arguments based on the central 
condensation and on the time scale of evolution vis-a-vis their 
observed numbers indicate that these BO. 5 - B2 giants are in the 
hydrogen-burning phase. Model sequences are constructed for 
stars of 15 and 20 to the end of hydrogen-burning. The pul- 
sational characteristics are then obtained by perturbing the 
stable models in the usual adiabatic approximation. Since the 
pulsational amplitude is greatest in the outer envelope, where 
the effective polytropic index is rapidly varying, it is shown 
that previous arguments based on polytropic models are liable 
to considerable error. Indeed, at some central hydrogen abundance 
which is higher for the lower masses, the periods and their 
ratios seem to be accounted for. Ideally, these quantities im- 
ply a unique mass and mean molecular weight for each observed 
star or} the H-R diagram, but comparison of theory and observa- 
tion gives at present merely a mass range of 10 - 20 M^ and 
probably a "normal" chemical composition. At any rate, the 
hypothesis of evolution of main sequence 09 - B1 stars across 
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the instabil-ity strip seems to be correct. / Uncertainties in 
the semiconvactive theory of massive stars would appear to be 
irrelevant, since the cpaasi-stable zone has almost no effect on 
the pulsational eigenfrequencies . Theoretical .and observational 
arguments indicate that the effects of rotation and mass loss will 
be small. The period-lximinosity relation, observed luminosity 
classes, and location of the Trumpler turn-off in clusters tend 
to confirm our prediction that the 0 Cephei stars of lower mass 
fall closer to the initial main sequence. The Wolf-Rayet objects 
form an apparent extension of the instability strip to higher 
masses (0 stars) . 
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I . INTRODUCTION 


The 0 Cephei stars are a group of short-period, pulsating 
variables. They form a well-defined sequence on the H-R diagram, 
running closely parallel to, although somewhat steeper than the 
main sequence at BO. 5 - B2 and having luminosity classes III and 
IV. They obey a period- luminosity law, with periods ranging from 
3.3 to -7 hours (toward higher luminosities). The observational 
data on the 0 Cephei stars, also called 0 Canis Majoris variables, 
have been summarized by Struve (1955j^) . There are now 18 or so 
known members of this group (van Hoof 1962^) . 

Since some of the 0 Cephei stars exhibit a beat phenomenon, 
caused by two nearly equal periods and appearing only in those 
members of the group with the highest projected rotational 
velocities (30 to 60 km/sec), Ledoux (1951, 1958) suggested that 
the instability was caused by rotation, and furthermore that it 
was probably manifested by non-radial oscillations. Basing their 
conclusions on crude interior models and only a rough notion of 
the mean densities of 0 Cephei stars, Ledoux and Walraven (1958) 
tried to show that the period-density relation calculated on the 
basis of radial oscillations was in conflict with the observations. 
More recently Gurm (1963) considered radial oscillations of a 
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main-sequence B1 star of 10 (Kushwaha 1957) and came to the 
same conclusion. Adopting the concept of non-radial oscilla- 
tions and a suggestion by Chandrasekhar and Lebovitz (1962) , 
Bohm-Vi tense (1963) showed that the degree of period commensur- 
ability^ could be explained by the observed rotational velocities 
and vice versa , if the 3 Cephei stars were observed pole-on. 

She used the latter condition to explain the relatively small 
number of g Cephei stars among early B giants. However, her 
conclusions were based on idealized polytropic models of rotating 
stars from Schwarzschild (1958) and only a rough application of 
results from the Chandrasekhar-Lebovitz theory of non-radial 
oscillations of rotating fluid spheres. 

Since rotation is not always observed in the g Cephei stars, 
and is in any case small in comparison with 100 to 200 km/sec for 
main-sequence and giant stars of the Scune spectral type (McNamara 
and Hansen 1961) , it seems worthwhile, still, to consider the 
variables as non-rotating stars undergoing radial oscillations. 
Struve (1955^) has emphasized that many of the g Cephei stars 
are relatively simple variable stars, completely analogous to 
other variables that are believed to be in purely radial oscilla- 
tion. In fact, van Hoof (1962^,) has shown that many of the 
observed features of the g Cephei phenomenon may be explained 
by the interference of several simultaneously excited modes of 



radial oscillation. ^ 

It is the purpose of this paper to determine at which 
evolutionary stage massive stars become 0 Cephei variables, and 
to see whether the relevant observations may be explained on, the 
assiomption of radial pulsation. 

II. PRELIMINARY EVOLUTIONARY CONSIDERATIONS 

From a discussion by Schmalberger (1960) , the 3 Cephei stars 
lie on the H-R diagram near the locus of models at the end of 
hydrogen -burning. This locus is defined by the first turn back 
toward the main sequence. Thus the 3 Cephei stars appear to 
have evolved from main sequence 09 - B1 stars. Since calculated 
luminosities for models of upper main sequence stars are fairly 
accurate, we may use the observed period-absolute magnitude re- 
lation to obtain the period-mass relation: II (hours) = 0.35 M/M^, 

where the mass range is roughly 10 - 20 M^. The major uncertainty 
lies in the bolometric corrections to the observed magnitudes. 

There are no observed masses of these stars. 

/ 

'If the 3 Cephei stars are indeed situated along the locus 
of secondary contraction, then some ambiguity arises because the 
evolutionary track for models of massive stars swings back quickly 



to the right after the brief turnback, forming an S-shaped curve. 
The luminosity becomes only slightly higher, but the internal 
structure is grossly changed from a core burning configuration 
to a contracting core with a surrounding, hydrogen-burning shell 

A 

(Sakashita, Ono, and Hayashi 1959; Hayashi and Cameron 1962; 
Stothers 1963 and 1964, hereafter called "Paper I" and "Paper II", 
respectively) . The contraction phase was suggested by Reddish 
(Discussion, Reddish and Sweet 1960) . 

Two lines of evidence point strongly to the former configura- 
tion for the g Cephei stars. The first line of evidence is based 
on the observed statistics of these stars. If, our sample of 18 
stars (van Hoof 1962jg.) is complete within a kiloparsec radius 
around the sun, in a Galaxy of effective R = 10 kiloparsecs, and 
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if there is a maximum of 1 x 10 0 - B2 stars in the Galaxy, 
then the number ratio of g Cephei stars to all 0 - B2 stars will 
be less than 1/50. In fact, McNamara and Hansen (1961) have ob- 
tained the ratio 1/9. Now according to Hayashi and Cameron (1962) 
the gravitational (core) contraction phase of a star of 15.6 
is 1/200 as long as the hydrogen-burning phase. If hydrogen- 
exhaus-fion is counted with the gravitational contraction phase, 
the ratio becomes 1/80. These ratios are far larger than the 
observed ones. 


The second line of evidence concerns the period ratios. 



according to the work of van Hoof and the results to be derived 
in this paper. During the hydrogen-burning phase, the calculated 
and observed period ratios agree at some characteristic central 
hydrogen abundance, dependent on the stellar mass. However, 
increasing central condensation leads to more and more discrepant 
values for these ratios. 

Hence we conclude that we must seek models for the 3 Cephei 
stars in the hydrogen-burning phase of evolution. 


III. STABLE MODELS 
Basic Physics 

The general structure assumed for models of massive stars 

has been outlined in Paper I. Here we adopt the same assumptions, 

notations, and equations as before. The adopted masses are 15 

and 20 M^, since at 10 electron scattering' is no longer a good 

approximation to the opacity throughout the star. The initial 

composition is again taken to be 
i 

X = 0.70, Y = 0.27, Z = 0.03, X__ = Z /2. (1) 

e e e G^O e 

The parameters in the formula for nuclear energy generation are. 


in the present case. 
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37 X 10“ 



- 1 1 4-. 1 


X. > 

>c \ 0 ^ 

y «- I S' 

It 

— 10^.6 , 

(2) 


b) Integration of Equilibrium Equations 

The construction of models proceeds in the manner outlined 
in Paper I, However, the fitting was accomplished completely 
automatically by the computer, as follows. We denote the en- 
velope eigenvalue by a. Thus for the initial model a= log C, 
and for the inhomogeneous models a = X (the composition exponent) . 
The integrations of the outer radiative Zone I and the assumed 
radiative intermediate Zone III are performed as in Paper I, 
and then smoothly continued into Zone IV, the convective core, 
when (n + 1) , = (n + 1) Care must be exercised in selecting 

the proper fitting point, since (n “ (n + the 

semiconvective zone, too. Once inside the core, the inward 
integrations are tested on the U-q plane at each point. If 
U > 3 or du/dq > 0, the current value of q (say, q^) is noted, 
and the' whole integration is begun again from the surface (for 
the initial model) or from q^ (for the inhomogeneous models) 
with a revised estimate 'of a based on the smallness of q^. When 
q^ is less than a suitably small prechosen value, an extrapolated 
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"value of 3 is obtained, 
c 

The integration method is now altered as follows. From the 
last trial integration a fitting point is selected just inside 
the core; the value of 3 = is noted. Now two inward integra- 

tions are performed for a and a + Aa, in to 3^. Here the values 
of U and V are noted. Similarly, two outward integrations from 

the center are performed for 3 and 8 + A3 . From the two in- 

c c c 

ward integrations we may form 6U/6a and 6V/6a, and from the two 

outward integrations 6 U /63 and 6 V /63 , all evaluated at 3^. We 

c c f 

obtain improved values of a and 3 ^ by using the increments Aa and 
3 ^ derived from a solution of the two simultaneous equations, 

- V-.r. = 15 - 'If; 

(3) 

This method is iterated until Aa and Ag^ are suitably small. 

When Zones I, III, and IV become properly fitted in U, V, 
and g, the semiconvective Zone II is integrated and fitted to 
Zone III' as shown in Paper I. 

' Results for Stable Models 

Table 1 contains the essential results for evolutionary 
sequences of six models obtained for stars of 15 and 20 
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Coinparison ir.ay be rr.ade with the analogous se-quer.ce calculatec! for 

30 K in Paner I. It shoulc be noted that in this paper the 
© 

adopted values of and are those of Allen (1963) , whereas 
in Papers I and II Chandrasekhar's (1939) values were used. In 
all comparisons with other work, v/e shall renormalize luminosities 
and radii to Allen's values whenever necessary. 

Apart from the increasing importance of the semiconvective 
zone from 15 to 30 the only other point worthy of special 

mention is that the initial decrease of central density for stars 
of intermediate and high mass becomes negligible or even vanishes 
for very massive stars, at some mass between 20 and 30 (cf. 
also Henyey, LeLevier, and Levee 1959) . It seems that the nuclear 
energy generation is not quite sufficient to expand the central 
regions against their slow gravitational contraction. This, would 
be due to the lower value of , which leads to a less steep 
gradient of temperature and hence of pressure. 

IV. PULSATING iMODELS 

a) Pu 1 sat ion Ecmation ■ 

The equation describing small, radial adiabatic pulsations 
of a spherically symmetric fluid may be written 
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n 




-i- 

or 
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c. o 

fc- V 


rr 


n Tr ■'‘^r 


■{ 


n p) 
I iH 


= 0, (4) 


where § is the radial displacement, ct/ 2~ = E ^ the frequency and 
n the period of oscillation. We have also the adiabatic exponent 


r. - a . ( 'c -l) 

M " + 12 ^ J 


(5) 


for a mixture of perfect gas and radiation (Chandrasekhar 1939) . 
Equation (4) is the equation derived by Ledoux (1939) if the 
relative amplitude 6r/r is introduced in place of § = 6r. Since 
the coefficients of the pulsation equation have singularities at 
the surface and center, we have the following boundary conditions 


sr= -r,F (2^ 

Sr " s. “ O • 


3r 


= 


fr- k") 

(r = 0) 


( 6 ) 


It will be convenient to evaluate an explicit expression 
for 3(r^P)/Sr. Since the ionization zones of hydrogen and helium 
are of negligible extent in massive stars, we may take y ~ 
througliout the whole star. Then we obtain from equation (5)' 




7T, 


- 2-^ 


S-7 




(7) 
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irind from the equation of state 

^ ^jLt^T'). ( 8 ) 

With the help of the homology invariants V = - d In P/d In r and 
n + 1 = d In p/d In T, we obtain the desired expression 

jfplf (r,r) = . <9) 


where 


, n-(3Vn.-3) 


B = 


Iillzlll 

yr, 


(10) 


Introducing equation (9) , V, and gP = kPT/tiH into equation 
(4) , we rewrite the pulsation equation as 




I 


This is the final form of the pulsation equation, taking radiation 
pressure exactly into account. In the case of no radiation 
pressure (1 - g = b = 0) , = Y and equation (11) reduces to 

the equation used by Stothers and Schwarzschild (1961) . 

In^- terms of the non-dimensional envelope variables defined in 
Paper I equation (11) becomes 



SI. 


2 . 







( 12 ) 
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wheye 

a z 

^ (13) 

2 . 

is the eigenvalue of the problem. For homologous stars (ju is 
•* 

constant, so that we obtain the familiar relation II /^ = const. 

. In the core of the star, the pulsation equation may be 
written as equation (12) with the envelope variables simply re~ 
placed by the starred core variables of Paper I. In this case 
the eigenvalue is 


u> 


MX S 





3(>c. 


(14) 


h) starting Series 

Because of the singularity at the stellar surface, the 

solution of the pulsation equation must be expanded in a power 

series around r = R. To do this, we first introduce the zero 

surface expansions from Paper I into equation (12) . We note 

the necessary expressions for the homology invariants: 

/ 

^ (15) 

« 

Then equation (12) becomes, near the surface. 
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^7< 


(16) 


where 


Ml 




r,‘ 


y= = 2-8 (2-y^) 


(17) 


and r^ = (g^) . We now expand f in a power series in terras of 

the small (negative) quantity x - 1: 


Jc = 0 ' 


Jt. 


(18) 


Introduction of this expression into equation (16) yields the 

2 

following recursion relation among the coefficients I (lu ): 

iC 

ioDc+s') =1 5 j^_3 + (3^; + 3JL-JJ-) - 

. > V. 

+ (t^s - r* ^ ^jb-i > 

where = |__2 “ §_2 ~ Note that the magnitudes of the co- 
efficients are relative, to the magnitude of Normalizing the 

amplitude at the surface (x = 1), we have that § =1. 

o 

The solution of the pulsation equation is similarly expanded 

in a power series around r = 0. From the center expansions of 

Paper I,-, we note that V is a quantity of second order and g « g , 

c 

so that 


XT = O j *V*= O ^ Y>-+ I = 


(2Q) 
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Then equation (12) becomes, near the center, 


1 ; 


— ~ 0 , 


( 21 ) 


which is independent of 3 and w . The solution is simply 




( 22 ) 


c,) Integration of Pulsation Equation 

It is convenient to split the second-order differential 

equation (12) into two coupled, first-order equations with the 

help of the auxiliary quantity '«}/ = S§/5x. By a suitable choice 

of independent variables in place of x for the various regions of 

the star, these two equations may be accurately integrated rather 

2 

close to the stellar center. Since the solution for arbitrary tu 

will, in general, diverge as the center is approached, the following 

procedure was used to obtain the correct eigenvalue that makes 

S and small near the center. For each mode, two trial values of 
2 

w are guessed, and integrations are performed in both cases down 

to a small, prechosen value of q (say • Then an improved value 
2 

of to is calculated from the magnitude and sign of 'll at q^. The 
procedure is iterated until i]; « 0 at q^. 

Improvement is now obtained by the same method that was used 
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in calculating the stable models (Section Illb) . We fit the 
core envelope solutions at using 5§/Sg and 3 §/3g as 

fitting parameters. The inward and outward values of these 
parameters are normalized by the ratio of the inward and outward 

values of § at An alternative procedure is to use equation 

* 2 2 
(14) for (ju directly from the trial value of uu and our knowledge 

of P and from the stable model. Then only one equation of (3) 

2 *2 

has to be used to obtain improved values of uu and hence of cu 

*2 2 

We note that, besides g^, lu depends on u) only through and 
t^/t* (from the stable model). 

^) Results for Pulsating Models 

The pulsational characteristics of massive stars are col- 

2 *2 

lected in Table 2. It is to be noted that, since uu and ou are 

independent of the nuclear energy generation and hence of the 

stellar radius, so are and the period ratios. They depend 

only on g^ and the total luminosity, apart from the age and 

initial chemical composition. The period-root-mean-density 

relatipn is given by Q = li/" (P/P^) • It does not change much from 

mass to mass, for models at the same evolutionary state. The 

same statement is crue -for P^/P and the period ratios. 

2 

We note that uu does not depend very strongly on' the central 
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condensation. For instance, between models 1 and 2 for 15 

? /P doubles. Likewise, between models 3 and 5 it doubles again, 
c 

2 . 

However, the change of oi xn the latter case is much smaller than 

2 

in the former. The reason is that cu depends on the rate of de- 
crease of S/ which drops rapidly in the outer envelope, and there- 
fore is insensitive to the precise interior conditions. Figure 1 
shows this in the case of the five calculated modes for model 4 

of 15 M . Since § is determined chiefly by the structure of the 

2 

envelope and hence by the total lumxnosxty, so xs uu . Now we 
note that the luminosity increases much less between models 3 and 
5 than between models 1 and 2. 

2 

Tables 1 and 2 also show that (ju and the period ratios de- 

o 

crease with p and F^. This may be seen by comparing the initial 
models for the three masses. The effect is similar to that found 
in the standard model (Ledoux and Walraven 1958) . However, de- 
tailed comparisons using the standard model appear to be very 
misleading, because the effective poly tropic index throughout the 
envelope actually drops far below n = 3 (see Figure 2) . This leads 
to grealtly different values of the pulsational quantities. In the 
last column of Table 3 , Schwarzschild ' s (1941) results on the 

standard model are listed. We see that our values of Q are even 

o 

larger than Q of the standard model, which Ledoux and Walraven 
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_found to be in excess of their "observed" value for the p Cephei 

stars, Q = 0.027. However, the values they used for the radii 
o 

of these stars were too large. 

Gunn (1963) calculated the pulsational eigenvalues for Kushwaha’s 

(1957) initial main sequence model of 10 and X =0.90, and found 

O e 

Fi^ = 2.12 hours. Since we see from Table 2 that the period more 
than doubles along the evolutionary track of massive stars, his con- 
clusion that our present theory of stars on the upper main sequence 
seems inadequate, because of the period disagreement, appears to be 
unwarranted. In fact, as we shall see in Section Via, the periods 
(but not their ratios) may be considerably altered through small 
changes in the chemical composition. Moreover, we have seen that 
arguments based on changing the central condensation of models do 
not lead very far, since Table 3 shows that even the standard model 
gives a fair representation of P /^. Finally, a gradient of ef- 

I 

fective polytropic index, such as occurs in the models calculated 

2 

here, seems necessary to yield correct values of p /p” and csj . 

c o 

Reddish and Sweet (1960) interpreted Struve's tentative sug- 
gestion of a secular period change in 3 Cephei in terms of the 
expanding radius during hydrogen-burning. Their rough result that 
the rate of increase is an order of magnitude smaller than that re- 
quired by observations -is confirmed by our detailed models. This 
is a further indication that beat phenomena ' and slow period changes 
in ; Cephei stars may be explained by interacting modes. 
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V. DATA ON THE BETA CEPHEI STARS 

I 

In Table 4 the results of van Hoof’s analyses of the light 
curves of five 0 Cephei stars are presented. A similar analysis 
of Canis Majoris is not presented, since it was essentially an 
interpolation (van Hoof 1963). However, these stars are repre- 
sentative of all the 3 Cephei stars, in that they include a broad 
range of periods, rotational velocities, and velocity amplitudes. 
The data of Table 4 were used to interpolate the models presented 
in Table 3. All the period rarios agreed; therefore we listed 
only the adopted each mass. 

To plot the observational data for 3 Cephei stars on the 
theoretical H-R diagram, we have used the list given by van Hoof 
(1962^) and the relations between spectrum, effective temperature, 
and bolometric correction given by Harris (1963). Harris's bolo- 
metric corrections were used for all luminosity classes (II - IV) . 
The luminosities were normalized by using 4.72 for the 

sun (Allen 1963). The resulting values of luminosity and effective 
temperature are not very different from those obtained by Schmal- 
berger (1960) . 

The evolutionary tracks on Figure 3 came from Papers I, II, 
and the present paper for 15 - 20 - 30 and from Kenyey 
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(1959) for 11 - 20 M^. Crosses denote the interpolated models 
of Table 3. 

It is clear from Figure 3, as well as Table 3, that the 
locus of constant period ratios is not the locus of secondary 
contraction. We shall not concern ourselves that the constant- 
ratio strip falls to the left of the observed g Cephei strip and 
that the periods do not quite agree; both of these discrepancies 
may easily be removed by a change in the radius, accomplishable 
in several ways, as shown in the next section. It is therefore of 
somewhat academic interest to recall that Struve (1955^) had 
earlier proposed evolution of stars op the g Cephei strip, in 
the absence of any detailed models at the time. 

VI. EFFECT OF CHANGING PHYSICAL PARAMETERS 


Chemical Comoosition 


In general, a change in the initial chemical composition of 
a star whose opacity is dominated by electron scattering will 
shift its evolutionary track along a diagonal line almost parallel 

to the initial main sequence in the H-R diagram. For, since 

4 * 1/2 

L ~ u and R ~ u / we have that T ~ p 
•^e ^e . e e 
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If we assume that all the 3 Cephei stars have the same mass 

and pulsate at the same evolutionary stage, but have different 

3/2 3/2 

initial chemical compositions, then n ~ R ~ p . Since U xs 

6 

observed to vary by a factor of 2, we must conclude that also 
varies by at least the same factor. This conclusion seems to be 
unjustified by observations of the upper main sequence. Moreover, 
the period ratios will be different at the same evolutionary stage, 
in contradiction of the observed rough constancy of these ratios. 
Finally, the variation in X^ necessary to produce the observed 
variation in luminosity is much more than a factor of 2. 

The last argument would also rule out the possibility of 
pulsation at differing Revolutionary stages, even though supple- 
mentary calculations indicate that the constant-ratio strip on 
the H-R diagram for stars of the same mass but differing initial 

chemical composition falls in the same way as in Figure 3. The 

2 

reason for the similar position of the strip is that and . the 

period ratios at the same evolutionary stage (1^) are smaller for 

stars of lower X (and hence higher luminosity) . 

e 

If, however, we assume with Schmalbesrger (1960) that all the 
0 Cephei stars fall along the locus of secondary contraction 
(X^ « 0.03) on the H-R- diagram, we can compute the necessary 
initial .chemical composition for each mass on the basis of constancy 
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of the period ratios. From Table 3 the model for 30 with 

X =0.70 almost fulfils the criteria. Using it as a standard, 

0 

we invoke homology arguments to calculate for other masses. 

From Paper I, as long as is not too low, the dimensionless 

structure of the star is specified only by the parameters A and 

C at a given evolutionary stage (1 ) , since the variable j may be 

c 

0 235 

replaced by 1 ' (Schwarzs chi Id 1958) . Then, holding A constant, 

we calculate that X must be 0.46 and 0.32, for 20 and 15 

respectively. Since the pulsational eigenvalues and hence period 

ratios depend only on the dimensionless structure of the star, 

these values of X^ seem to be necessary for the hypothesis of 

secondary contraction. They are unrealistically low, and should 

be lower still because the mass at which X^ = 0.70 is actually 

greater than 30 M^, and because the decreased values of X^ imply 

an increased luminosity. Therefore, to give agreement with the 

observed luminosities, even smaller masses and hence lower X would 

e 

have to be taken. Moreover, there is no reason to believe that the 
calculated periods would agree with the observations. In any case, 
it is .difficult to see why 3 Cephei stars of lower mass should have 
lower initial hydrogen abundances. 

Ideally, if the 3* Cephei strip were sufficiently well defined 
observationally , the stellar masses and chemical compositions could 
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be determined with greater accuracy. For a given model, specified 

by M and a , the evolutionary track crosses the strip on the H-R 
e 

diagram at a certain point, where both and the period ratios 
must agree with the observations. Since depends essentially 
on the stellar radius and the period ratios on the luminosity, the 
required model is therefore uniquely determined. 

b) Nuclear Enemy Generation 

The radius of massive stars is essentially determined by the 

3/2 

rate of nuclear energy generation. Since II ~ R , it is clear 
that we may seek agreement with the observed periods by adjusting 
or 3 .S well as by changing (see previous subsection) . 

To change II of the model for 15 M. that best fits the observa- 

, O 0 

tions of period ratios, from 3.7 to 5 days, we require an increase 
in e or X by a factor of 45. Such an increase seems inadmis- 
sible. Moreover, the increase does not scale homologously with 
inverse mass but is constant, so that agreement with observations 
at 15 produces too great a at 20 M^. It should also be re- 
called that changing R and hence the periods does not change the 


period ratios. 
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£.) Onacitv 

The inability of a reasonable change in the nuclear energy 
generation rate to produce agreement with the observed periods 
will not be disastrous, however. In our models we neglected 
opacity sources other than electron scattering, and it is certain 
that bound-free absorption processes will contribute non-negligibly 
in the outer envelope. The model sequence for 20 with = 0.68, 
computed by Henyey, LeLevier, and Levee (1959) , included these 
processes, and the resulting evolutionary sequence lies on the 
H-R diagram parallel to our sequence at very nearly the same lumi- 
nosity. It is, however, displaced to lower effective temperatures 
by an amount equivalent to a change in log (R/R^) equal to 0.06, 

after allowance for differences in X , X„^, and e . This change 

e CNO o 

is brought about almost directly by the opacity, since we have 
that R ~ H, from a dimensional analysis of the equation of energy 
transport by radiation. The change produces an increase of 11^ from 
5.2 to 6.4 days, more closely in agreement with observations. 

The question arises whether inclusion of bound-free absorption 
will change the pulsational eigenvalues. Undoubtedly it will to 
some extent, but the ratios of the modes, especially those of the 
higher modes, should remain fairly constant, because they are nearly 
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independent of the radius. We recall that it is these ratios 
that essentially fix the 3 Cephei striy . the evolutionary H-R 
diagram. 

It, may easily be shown, however, that the inclusion of bound- 
free absorption, which has an increasing effect at lower masses, 
actually serves to offset the line of constant period ratios 
farther from the locus of secondary contraction. Its inclusion is 
roughly equivalent to increasing X in the electron scattering opacity. 
As discussed in Section Via, the period ratios then also increase. 
Therefore, at lower masses, the model for which the ratios agree 
with observations lies closer to the initial main sequence. How- 
ever, as stated above, 'the effect should be small. 

,d ) Semi convection 

Since the extent of semiconvection in a star depends mainly 

on the luminosity, a lowering of the initial hydrogen abundance in 

stars of a given mass increases the eimount of semiconvection, 

4 

through the relation L ~ However, the semiconvective zone, 

as we ,have treated it, is unable to alter the pulsational eigen- 
values to a perceptible degree, even in the last hydrogen-burning 
model of a star of 30 M^. Hence uncertainties in the semiconvec- 
tive theory will probably not be reflected in the pulsational 
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characteristics of massive stars. For determination of these 
characteristics, it is adequate merely to consider the intermediate 
zones as wholly radiative. 


VII. DYNAMICAL EFFECTS 
Rotation 

Rotational. effects should be small, since McNamara and Hansen 
(1961) have shown that the average equatorial velocity for 3 Cephei 
stars is only 22 km/sec. In addition, theoretical evidence exists 
that large rotation may inhibit pulsations (Steinitz 1964) . How- 
ever, the fact that somewhat higher rotational velocities are ob- 
served in the case of 3 Cephei stars exhibiting a beat phenomenon 
indicates that the presence of secondary periods may be due to the 
rotation. At any rate, the work of McNamara and Hansen (1961) and 
the summary by Struve (1955^) strengthen the case for purely radial 
oscillations in the slowly rotating 3 Cephei stars. In view of 
the observed period-luminosity relation and van Hoof's work, 'this 
statement probably holds true also for the more rapidly rotating 


stars. 
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Energizing Mechanisms 

Schwarzschild and Harm (1959) have shown that upper main 
sequence stars become pulsationally unstable above about SO M^. 

At these masses approaches 4/3 because of the high radiation 
pressure, and nuclear reactions are able to supply enough energy 
to maintain the instability. Stars of lower mass should be pulsa- 
tionally stable nuclear-wise, and explicit calculations by Gurm 
(1963) for 10 confirm this expectation. 

, Rotational energy is a second possible source for the pulsa- 
tions. As evolution off the main sequence proceeds, the envelope 
expands and hence displays a lower rotational velocity. At the 
same time, the core is - contracting and, in the absence of dis- 
sipative mechanisms, will rotate faster. Since the 8 Cephei stars 
as a group show rotational velocities far less than the velocities 
shown by stable stars of the same spectral and luminosity classes, 
possibly the angular momentum has been transferred to the core and 
used in energizing the pulsations. However, Figure 1 suggests that 
rotational energy of the core may not be a source for the pulsations 
since the pulsational amplitude here is less than 1 per cent of its- 
value at the surface. It is partly for the same reason that nuclear 
energy fails as a source. 

Further, as we have already mentioned, indirect theoretical 
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evidence shows that rotation of a star as a whole tends to in-, 
hibit pulsations (Steinitz 1964). In analogy with the cepheids, 
if the hydrogen or helium ionization zones are the source of the 
pulsations, then they originate in the atmosphere or not far be- 
low the^ photosphere. Atmospheric expulsion is indeed what is re- 
quired to explain the beat phenomenon according to the theory of 
Struve and Odgers (Struve 1955]^) . If molecular dissociation is 
the cause of the instability, then some reasonably abundant molecule 
must be dissociating ate temperatures typical of BO. 5 - B2 giants. 

I 

i £_) Mass Loss 

I 

If the theory of Struve and Odgers (Struve 1955^) is correct, 
the 3 Cephei phenomenon may be explained by the ejection, decelera- 
tion, and subsequent inf all of an atmosphere. In any case, it is 
to be expected that soii^fe mass will be lost (Sahade in Discussion, 
Reddish and Sweet I960)'. We should like now to examine whether the 
3 Cephei strip (or constant-ratio strip) is actually the evolu- 
tionary track of a star losing mass. ; 

If the star remains chemically inhomogeneous, a constant- 

ratio line cannot be maintained since X along the line increases 

c 

as the stellar mass is [lower (Table 3). If, however, the incipient 
instability causes and (then maintains complete mixing of the 
stellar material, a corj(Stant-ratio line' might be maintained since 


30 - 


(= X^) decreases with the mass. For average values of M and L 
taken from Table 1, the lifetime of ' a star of initially 20 to 

7 

reach 10 will be At = E AX <M>/<L> = 10 years. Hence the mean 

—6 

rate of mass loss is 10 M^/year; this rate might not be unreason- 
able. 

Three arguments seem to rule out complete mixing, however. 

-2 

First, the mass would be forced to decrease with X as u in order 

e e 

to preserve constancy of the period ratios (see Section via) . 

4 3 

Second, since L ~ u M|/ (1 + X ) ~ M/ (1 + X ) , and both M and 

e e e 

1 + X^ decrease by about the same factor, L will not change very 
much, in contradiction to the observations. Third, complete mixing 
not only restores stars to the initial main sequence, but as hydro- 
gen is consumed, it produces a track to the left . A compromise 
based on partial mixing may be ruled out by the same argument ap- 
plied against the inhomogeneous case. 

We conclude that theory predicts little mass loss, and in the 
absence of any direct observational evidence to the contrary, we 
have assumed that stars must therefore evolve across the instabi- 
lity s^trip. Since the strip is so narrow, the time scale of 
evolution across it must be small, and therefore the mass loss in 
any case will be small r 

■1 ' 
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VIII. INTERPRETATION OF THE H-R DIAGRAM 

We should now like to see whether the rapid drop in lumino- 
sity along the constant-ratio line is more compatible with the 
observations than the gentler drop occurring strictly along the 
locus of secondary contraction, as suggested by Schmalberger (1960) . 
First, since no known 0 Cephei stars are members of a binary system, 
we have relied on the model calculations to place the mass limits 
at 10 and 20 M^, roughly. Then our constant-ratio models predict 
a period— luminosity law II ~ L * in this range. (Extrapolating 
from 15 to 10 M^, we should actually have an exponent slightly 
less than 0.40.) Since II will change by a roughly constant multi- 
ple for all model masses if the constants determining R are changed, 
the exponent 0.40 will! remain unchanged for horizontal shifts of 

the evolutionary tracks in the H— R diagram. Now all the 0 Cephei 

0. 25 

stars taken together (van Hoof 1962^) yield a law II ~ L . How- 
ever, only four of them have accurately determined luminosities . 
These are members of the Scorpio-Centaurus cluster, and include 
0 Ophiuchus and 0 Crucis from Table 4. They yield the law 
H ~ L . Good agreement is therefore found with the theoretical 
law. r 

Second, the luminosity class drops from II - III for the 


M 
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variables of earliest spectral type (BO. 5) to IV for those of 
latest spectral type (B2) . This suggests that the 3 Cephei strip 
does indeed approach the main sequence closer than does the locus 
of secondary contraction, which should probably not show a drop, 
or at least a large one, in luminosity class. We note further 
that the magnitude difference between class IV and V stars attains 
a minimum at B2 (Arp 1958) . From Figure 3 the constant-ratio 
strip, extrapolated, would run close to the initial main sequence 
at 10 (B2) . 

Third, observations of early-type clusters and associations 
indicate that the tip of the Trumpler turn-off, which by age argu- 
ments is believed to represent the point of secondary contraction, 
occurs at luminosity class III. For example, in I Geminorum the 
turn-off from the initial main sequence appears at Bl V and the 
tip of the turn-off at Bl III (Crawford, Limber, Mendoza, Schulte, 
Steinman, and Swihart 1955) . This suggests that the Bl IV 3 Cephei 
stars would not have reached the end of hydrogen -burning. 

The observed spectral (or mass) range of the 3 Cephei stars 
is remarkably well defined. One indication that the lower mass 
cutoff should occur near B2 is that our calculations show an in- 
tersecuion of the extrapolated constant-ratio strip with the main 
sequence near 10 M^. Secondly, noting that the 3 Cephei stars 
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occur only among sharp-line (slowly rotating) early B stars, 

I 

McNamara and Kansen (1961) ascribe the cutoff to increasing rota- 
tional velocities among the late-type B stars. (This increase is 
observed in both luminosity classes V and III [Allen 1963].) Thus 
observations, as well as theory, indicate that large rotation tends 
to inhibit pulsational instability. 

An unfruitful suggestion regarding the uoner limit to the 
mass is that semiconvection starts to become important in stellar 
envelopes at about 20 M^. Although convection tends to damp pulsa- 
tions, the semiconvective zone is too ineffective and lies too deep 
for this purpose (see Tig. 1) . 

Observational evidence exists, however, for the continuation 
of instability up to the highest masses. But in the case of the 
O giants, the instability manifests itself in the Wolf-Rayet 
phenomenon. Westerlund (1961) and Westerlund and Smith (1963) 
have shown that in the-H-R diagrams of 0 clusters in the Large 
Magellanic Cloud, the Wolf-Rayer stars appear invariably at the 
tip of the Trumpler turn-off. These authors suggest that the masses 
of rhe; Wolf-Rayet starp lie between 20 and 60 which is what we 
require to explain them as an "extension" of th.e 3 Cephei strip, 
now occurring close torthe locus of secondary contraction. 

Sahade (1962) gives a table of computed masses for some of 
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the galactic Wolf-Rayet stars. Although they appear to be less 
massive than their OB companions, the three luminosity classes 
given for the companions are all class I. Hence we expect these 
companions to be more massive, since they have presumably evolved 
further (past the Wolf-Rayet phase of pseudo-class 0 III) . 

Finally, in an analysis of Wolf-Rayet spectra Smith (1955) 
has reported variability of emission line intensities on a time 
scale of a few hours, in analogy with the Of stars (Oke 1954) . We 
emphasize, however, that the observable form of the instability and 
probably the energizing mechanisms sustaining it are wholly dif- 
ferent for the Wolf-Rayet and 3 Cephei stars. Why an apparent 
changeover should occur at BO is unknown. 

IX. CONCLUSION 


It appears that radial pulsations of hydrogen-burning giants may 
be adequate to explain the observations of p Cephei stars. Non- 
radial or strictly atmospheric pulsations need not be invoked. How- 
ever, ;the effect of rotation and possible mass loss cannot be defi- 
nitely ascertained, although it is probable that they are small. 

In any case, theory cannot yet say why the p Cephei stars lie in 

i 

the distinct range BO. 5 - B2 III - IV, nor why only some stars in 
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this range become variable, nor what the sources maintaining 
several simultaneously excited modes may be. 


A preliminary theoretical discussion of the 3 Cephei stars ap- 
peared in the author's doctoral dissertation, Harvard University 
(1963), which was supported in part by a Harvard scholarship during 
the first half of the academic year 1963 - 1964. Another part of 
the work reported in this paper was supported by an NAS-NRC Post- 
doctoral Resident Research Associateship under the National Aero- 
nautics and Space Administration. It is a pleasure to thank 
Dr. Leon Lucy for discussions and Dr. Robert Jastrow for his 
hospitality at the Institute for Space Studies. 
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Table 3 


Interpolated Model Characteristics for the Phase 
of Equal Calculated and Observed Period Ratios 



' 15 M 

O 

20 

© 

30 Mg 
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log (L/L^) 

4.57 

4.94 

5.40 


log T 

e 
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4.55 

4.56 

1 

^4 

0.33 
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^ 6 

T (10 years) 

5.6 

5.2 

4.6 

.... 

log{R/R^) 

0.75 

0.89. 

1.10 

* • • • 
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X 
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3.7 

5.2 

9 


Qq (days) 
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0. 046 
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0.038 

nyn 

1 o 

0.675 

0.678 

0.68 

0.738 



Table 4 


Observed Period Ratios in 3 Cephei Stars 



9 Oph 

V Eri 
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3.37 
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4.57 
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. 
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nVn 
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• 
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II-/II 
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0.336 
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.-13/n^ 
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0. 599 
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0.602: 

0. 599 

n/Hi 

0.501: 
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0.499 

0. 501: 

0.496 

van Hoof 

1962b 

1961 

1962c 
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FIGURE CAPTIOUS 


Fig. 1 


Fig. 2 


Fig. 3 


- Normalized pulsation amplitude as a function of radius 
fraction for model 4 of 15 M^. Solutions are labeled 
with the mode number. Roman numerals designate the 
stellar zones (Paper I), which are marhed off by 
vertical lines. 

- Effective polytropic index as a function of mass 
fraction for the standard model and for models O and 
4 of 15 M . A dot marks the boundary of the convective 
core. 

. - Theoretical H-R diagram of the upper main sequence, 

including observational points for the 3 Cephei stars. 
Model sequences for 15-20-30 M^ are cue to Stothers 
and those for 11 - 20 to Kenyey, LeLevier, and Levee 
Crosses mark the interpolated models for equal cal- 


culated and observed period ratios. 
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